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ABSTRACT 


An incompressible boundary-layer stability analysis of Laminar Flow 
Control (LFC) experimental data was completed and the results are presented. 
This analysis was undertaken for three reasons: to study laminar boundary- 

layer stability on a modern swept LFC airfoil; to calculate incompressible 
design limits of linear stability theory as applied to a modern airfoil at 
high subsonic speeds; and to verify the use of linear stability theory as a 
design tool. The experimental data were taken from the slotted LFC experiment 
recently completed in the NASA Langley 8-Ft. Transonic Pressure Tunnel. 

Linear stability theory was applied and the results were compared with 
transition data to arrive at correlated n-f actors. Results of the analysis 
showed that for the configuration and cases studied, Tollmien-Schlichting (TS) 
amplification was the dominating disturbance influencing transition. For 
these cases, incompressible linear stability theory correlated with an 
n- factor for TS waves of approximately 10 at transition. The n- factor method 
correlated rather consistently to this value despite a number of non-ideal 
conditions which indicates the method is useful as a design tool for advanced 
laminar flow airfoils. 
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INTRODUCTION 


Renewed interest in improving airplane performance characteristics has 
forced researchers to re-examine the possibility of achieving extensive 
regions of laminar flow on modern airplane wings. Aerodynamicists have long 
known that reduction in drag and an increase in performance would be realized 
if laminar flow were attained during cruising flight. The performance and 
feasibility of both natural laminar flow (NLF) and controlled laminar flow 
(LFC) airfoils were tested in the forties and sixties with mixed results. 
Though for some of these cases full -chord laminar flow was achieved, it could 
not be maintained at flight conditions representative of current modern 
transport designs. The construction techniques of the day did not produce a 
sufficiently smooth wing for laminar flow at flight Reynolds numbers. 
Furthermore, the poor off-design flight and stall characteristics were 
unimpressive. These results led to the dismissal of laminar flow as a viable 
drag reduction technique. However, recent advances in technology, along with 
a marked increase in fuel costs have sparked renewed research in laminar flow 
airfoils. New design methods as well as new construction materials and 
techniques have increased the viability of such concepts. 

The design of laminar flow airfoils is dependent upon tools to model the 
boundary-layer characteristics correctly. Accurate prediction of the 
transition location is essential for determining the extent of laminar flow. 
Currently, the most reliable prediction technique is through the use of the e n 
method, which is based on linear stability theory. The e 11 , or n-factor, 
method is a semi -empirical technique which has gained credence from the fact 
that transition tends to correlate rather consistently with an n-factor value 
depending on the disturbance environment. For low disturbance tunnels, 
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n-f actors of n i = 9 corresponds fairly well to transition as established by 
Smith (ref. 1) and Van Ingen (ref. 2). For free flight, where turbulence 
levels are lower, transition seems to occur at much higher n-f actors . 3,4 
However, at present time there is only a limited amount of data available at 
flight conditions which are typical of high-performance transports. 

Validation of the use of the n-f actor method for the higher Mach and Reynolds 
numbers of typical transport flight conditions is needed. 

An experiment designed to examine supercritical LFC technology on a swept 
wing was conducted by NASA personnel in the Langley 8-Ft. Transonic Pressure 
Tunnel. 5 The overall objective of the experiment was to evaluate and 
document combined suction laminarization and supercritical airfoil technology 
at conditions which are typical of high-performance transports. Since linear 
stability theory was used in the design of this airfoil, an underlying goal of 
the experiment was to validate the use of this method as a design tool in the 
transonic regime. Therefore, the current thesis research, a stability 
analysis of the LFC experimental data, was conducted to correlate n-factors 
with transition at these high-performance transport flight conditions. The 
amplification of small disturbances for the swept LFC airfoil was calculated 
using the experimental pressure and suction data. Experimentally determined 
transition locations were compared with the results of these calcualtions and 
allowable amplification ratio limits were established. Measurements of the 
tunnel noise levels were made in order to assess the disturbance 
environment. The purpose of this thesis is to examine the disturbance growth 
problem in order to determine the design limitations of incompressible linear 
stability theory. 
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LINEAR STABILITY THEORY 


Historical Development 

The accurate prediction of transition is an important and formidable 
problem facing researchers today. A physically correct model of the tran- 
sition phenomenon is, as yet, non-existent because of the complex nature of 
the problem. Currently, the most reliable method available is based on linear 
stability theory. Although this method has met with much success, it has 
received some criticism because of the limitations of the model itself. A 
discussion of these limitations follows in a later section. A brief history 
of the development of linear stability theory and its application as a 
transition prediction tool is given next. While a large number of scientists, 
mathematicians, and researchers have played an important role in the 
development, only a few of the more noted contributions will be discussed. 

This section is intended only to outline the more important steps in the 
development of a reliable transition predictor. 

Much of the groundwork for the development of linear stability theory was 
laid down at or before the turn of the century. This early work on boundary- 
layer stability provided a strong foundation with which to develop the current 
transition prediction technique. The first examinations of consequence into 
the transition process were made separately by Helmholtz (1868), Kelvin 
(1871), and Rayleigh (1880). Their examination of the instability of inviscid 
flows provided some insights into the mechanisms of the transition phenomenon 
and also provided the initial theoretical formulation of the problem. Later, 
Reynolds (1883) performed his classic experiments on transition in a pipe. 
These experiments, which demonstrated the qualitative differences between 
laminar and turbulent flows, 6 provided the initial experimental investigation 
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into the problem of transition. An important breakthrough was made by Prandtl 
(1904) with his work on boundary-layer theory which allowed unification of 
experimental and theoretical work. Boundary-layer theory allowed experimental 
verification of the effects of viscosity. Prandtl (1921) also demonstrated 
the destabilizing effect of viscosity through the use of a simple linear 
model. 7 Previously, viscosity was thought to have a stabilizing influence. 

The next important step in the development of linear stability theory was 
work done on the stability equation and two-dimensional disturbance waves. 

The theoretical understanding of the stability problem took a decisive turn 
when Orr (1907) and Sommerfeld (1908) separately developed the stability 
equation which became known as the Orr-Sommerfel d equation. Although they 
succeeded in setting up the equation, they were unable to generate the 
solution to the problem because of its complexity. It wasn't until various 
heuristic methods of approximation, 8 and later adequately developed asymptotic 
analysis methods, were used that solutions to the Orr-Sommerfel d equation were 
found. Work by Tollmien (1929), Schlichting (1935) and others was 
instrumental in this regard. Their work resulted in a highly-developed theory 
for the stability of the (two-dimensional) Blasius boundary layer. 7 Tollmien 
and Schlichting's model of transition remained experimentally unverified until 
the brilliant experimental work of Schubauer and Skramstead (1948) 
demonstrated the existence of instability waves in a boundary layer. Their 
work provided conclusive evidence that instability waves exist in a boundary 
layer with small external disturbances, and that Toll mi en-Schlich ting (TS) 
theory accounted for this behavior. More importantly, Schubauer and 
Skramstead showed that instability waves were a necessary precursor of 
transition. 7 These two-dimensional instability waves were found to travel 
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roughly in the direction of the free-stream flow and were subsequently 
referred to as Tollmien-Schlichting (TS ) waves. 

The need for wing sweep required extension of linear stability theory 
into three dimensions. It was Gray (1952) who first noticed that swept-wing 
boundary layers became turbulent much closer to the leading edge than 
straight-wing boundary layers. 9 Earlier, theoretical work in three 
dimensions had been simplified by Squire (1933) with a transformation that 
reduced the three-dimensional problem to an equivalent two-dimensional 
problem. This meant that two-dimensional solutions could be applied to three- 
dimensional flows. Experimental investigation of three-dimensional flows by 
Gregory, Stuart and Walker (1955) found in the boundary layer the existence of 
a set of vortices with their axes lying roughly in the direction of the local 
external stream. These vortices are referred to as cross-flow (CF) 
vortices. This confirmed the existence of a different mode or mechanism of 
transition when dealing with the three-dimensionality of swept wings. Flight 
tests were also important in exposing the mechanisms of boundary- layer insta- 
bility due to swept wings. Pfenninger's (1958) work on swept wings with 
suction helped to provide an understanding of the basic problems of swept LFC 
boundary layers. 

The underlying goal for the development of linear stability theory was to 
develop a consistent and reliable method for transition prediction. Previous 
methods for predicting transition proved unsuccessful until Smith (1956) and 
Van Ingen (1956) independently developed a semi -empirical method based on 
stability theory. 10 This, so-called, e n or n-factor method has had a great 
deal of success in predicting transition and is currently the best known 
method of modeling the extent of laminar flow on an airfoil. Srokowski and 
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Orzag (1977) developed a user-oriented n-f actor computer program called SALLY 
which stands for Stability Analysis, Local, Linear, ^Incompressible. This 
program solves the temporal 0-S equation in order to find the amplification of 
disturbances within a boundary layer. Although this method ignores the 
influence of the initial disturbance environment, it has proved to be useful 
when calibrated against experimental data. A simplified code based on a table 
look-up method of finding the amplification of crossflow type disturbances was 
developed by Dagenhart (1981). Both codes have been shown to be adequate for 
indicating correct trends of the influence of various mean-profile modifiers 
upon transition, 11 * 13 thereby making them excellent design tools. 

n-F actor Method 

The n-factor method, which is based on linear stability theory has proven 
to be the most effective way of predicting transition at present. There are, 
of course, limitations to this method which are due to the assumptions imple- 
mented within linear-stability theory. An understanding of these limitations, 
and their possible implications on the results, is essential in order to 
properly apply the method and to get meaningful results. Therefore, a discus- 
sion of the relationship of linear stability theory and the n-factor method 
would be beneficial. First, a brief outline of the derivation of the 
stability equation, with emphasis on the important assumptions, will be 
presented. Then, discussion of the solutions of the stability equation as 
used in the n-factor method will follow. Finally, the assumptions of the 
stability theory will be examined in more detail in the hope of revealing its 
limitations. A more detailed explanation of the derivation of the stability 
equations can be found in texts by Betchov and Criminal e (ref. 12), Rosenhead 
(ref. 9), and/or Dagenhart (ref. 13). 
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The derivation of the stability equation to be used in this analysis 
begins with the incompressible, three-dimensional, time-dependent Navier- 
Stokes equations. A perturbation assumption is first used to simplify these 
equations. The perturbation assumption reduces the flow within the boundary 
layer to a superposition of mean-flow and perturbation components. The mean 
flow term is assumed to be steady, parallel, and laminar. The resulting 
perturbation equations with the mean removed are then linearized by neglecting 
the quadratic terms in the perturbation quantities. Separating variables by 
assuming a complex form of the perturbation functions reduces these equations 
to ordinary differential form. This representation models the perturbations 
as travelling periodic disturbance waves that have a certain frequency and 
wavelength. These disturbances are allowed to amplify or damp in either space 
or time. The general form of perturbation function is 

q(x, y, z, t) . +(y) e i(ax + 62 - “*> (1) 

wh ere q(x, y, z, t) represents any of the perturbation quantities. A further 
simplification is introduced through the Squire transformation, which reduces 
the three-dimensional disturbance equation to a single quasi -two-dimensional 
disturbance equation in the direction of propagation of the wave front. 13 
These transformations are 


a x = ax + ez (2a) 

a U = aU + 0W (2b) 


and 
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(2c) 


a 2 = a 2 + ft 2 . 


The angle between the direction of propagation of the disturbance wave front 
and the local inviscid flow at the edge of the boundary layer is defined as 
the wave angle. The mean velocity perpendicular to the wall is then neglected 
to obtain the Orr-Sommerfeld equation. This equation is the governing 
stability equation and can be used to model a two-dimensional travelling wave 
disturbance in a three-dimensional flow field. The Orr-Sommerfeld equation, 
given in reference 13, is listed here in its dimensionless form: 


(0 - ~) U" - ct 2 <j>) - tJ"4> = - (i/aR) ( <J>"" - 2$V + a^f) (3) 

a 


where <j> 

U 

O' 

U) 

a 

R 


is the disturbance amplitude as a function of the normal 
coordi nate 

is the steady-state velocity component in direction normal 
to disturbance wave front 

denotes differentation with respect to the normal coordinate 
is circular frequency of oscillation 

is the wave number corresponding to the direction normal to 
the wave front and along the wing surface 

is the Reynolds number. 


It can be shown that the amplitude of the disturbance q(x, y, t) grows 
spatially as 


A ~ 


r*. 



X 


(4) 


or 


/V /N. 



(5) 
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The growth rate a is defined as 


_ 1 dA _ - 

a = -j- — - - a- 

A dx 1 


( 6 ) 


If the local growth rate is integrated from the initial point of 
amplification, total amplification can be found as 


/* a dx = /* -j- ~ dx = nr) — = n. (7) 

x„ x„ dx o 

o o 

This defines the n- factor where n is the logarithm of the disturbance 
amplification from the neutral point. Using (6) in (7) yields 

n(x) = - J* a. (x) dx. (8) 

x 

x o 

Solutions to the stability (Orr-Sommerfeld) equation provide the 
necessary information for analyzing the growth of disturbances within the 
boundary layer. The Orr-Sommerfeld equation, along with the appropriate 
boundary conditions, constitutes a complex eigenvalue problem whose solution 
depends on the local mean flow field and the frequency, wavelength, and/or 
wave angle of interest. Two of the three variables (frequency, wavelength, or 
wave angle) may be specified with the third given by the solution of the 0-S 
equation. For these given conditions, a local disturbance growth rate can be 
found. Integrating the local growth rate, as in eq. 8, from the initial point 
of amplification (neutral point) gives the logarithm of the total 
amplification. The total amplification within the boundary layer is used as 
an indicator of transition. This is the basis of Smith and Van Ingen's n- 
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factor method where n is a running integral of the local amplification rate 
from the neutral point. The n-f actor method is based on the fact that at a 
given value of n, transition is likely to occur. Previous studies have 
indicated that correlation of the n-factor method with low-disturbance wind 
tunnel transition data occurs around a value on the order of ten. 1,2,10,11 

Several assumptions were used in the derivation of the Orr-Sommerfel d 
equation: (1) flow is incompressible, (2) flow is parallel, (3) nonlinear 

terms can be neglected, and (4) mean flow perpendicular to the wall can be 
neglected. A discussion of the effect of each assumption on the solution 
follows. The first assumption, that of incompressible flow, means that the 
density fluctuations are neglected. This, of course, is appropriate at low 
Mach numbers where the effect of compressibility is very small, but at the 
Mach numbers of the present analysis, compressibility effects should be 
considered. However, studies by Mack (ref. 19) and others have shown that the 
use of incompressible codes overpredict the local amplification rates for Mach 
numbers in the present range. Thus, compressibility has a stabilizing 
influence so that the use of an incompressible stability code should 
overpredict the disturbance amplification (a "conservative" approach). The 
second assumption used is that the boundary-layer mean flow is parallel. 
According to Schlichting (ref. 14) parallel flow is a good approximation of 
boundary-layer flow since the dependence of the velocity u on the y-coordinate 
is at least an order of magnitude greater than that on x. For the present 
analysis, the boundary layer is altered by suction through closely spaced 
slots. Distributed, or continuous, suction tends to make the boundary layer 
more nearly parallel after the initial growth in thickness near the leading 
edge. On the other hand, discrete suction (through slots) could cause local 
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"dimpling" of the streamlines which is counter to the parallel flow 
assumption. However, the effect of this streamline distortion should be 
small. Thus the parallel flow assumption is fairly good. The next 
assumption, linearization, means that the quadratic terms within the governing 
equations are neglected. Strictly speaking, this is justified only for 
infinitesmally small disturbances. In reality, the perturbation amplitudes 
are not infinitesmally small, but small enough such that the amplification 
process can be considered essentially linear. It is assumed that the linear 
amplification of disturbances dominates the process leading to transition and 
that nonlinear interaction and breakdown is reserved for a relatively short 
period just prior to transition. Thus, correlation of the n-factor method 
with transition data should be made at the beginning of transition (i.e., 
first appearance of turbulent bursts). The fourth assumption, that the normal 
velocity component is zero, seems contradictory to the fact that there is 
suction. However, as noted by Dagenhart (ref. 13), the most important effect 
of suction is modification of the streamwise mean flow velocity profile and 
its derivatives. Therefore, the impact of suction is manifested through the 
streamwise velocity terms retained in the Orr-Sommerfeld equation, and the 
error in neglecting the normal velocity component should be small. 

The theory as used in this research has been shown to be in reasonable 
agreement with experimental data. 13 There are discrepancies between theory 
and experiment, however, which can be attributed to these various 
assumptions. Independent work, such as Malik's work on compressibility and 
Saric and Nayfeh's work on non-parallel effects, is continuing in order to 
reintroduce these assumptions into the theory. Criticisms of linear stability 
theory, however, arises not from the use of these assumptions but from the 
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inability of the general theory to properly model the complete transition 
process. The mechanisms of transition can be thought of as occuring in three 
phases: (1) introduction of small disturbances into the boundary layer; (2) 

subsequent amplification of these disturbances by linear processes within the 
boundary layer; and (3) nonlinear breakdown to turbulence. The n-factor 
method and linear stability theory do a good job of modeling the second phase 
of transition, the linear amplification of small disturbances. The third 
phase, the nonlinear breakdown, is, of course, beyond the capabilities of a 
linear theory. However, if the beginning of transition is used as the 
correlation point, and if the interaction of streamwise and crossflow type 
disturbances are kept to a minimum, then this should not be important. The 
first phase, the introduction of small disturbances into the boundary layer 
(Morkovin's receptivity problem), on the other hand, is unaccounted for by 
linear stability theory. The influence of the external disturbance 
environment on boundary-layer stability remains an unknown. 

The disturbance environment appears to play an important role in 
determining transition n-factors. Studies by Pfenninger, Bacon, and Carlson 
(refs. 28-31) and Spangler and Wells (ref. 24) have shown the importance of 
the effect of acoustical noise on the stability of laminar boundary layers. 
These reports suggest that the spectral energy distribution near the critical 
band of amplified frequencies, and not the overall energy levels, is the more 
important influence upon transition. The relationship between transition and 
free-stream disturbance spectra, however, is not fully known at this time. 
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Boundary-Layer Instabi 1 i ti es 

Boundary-layer transition is a fluid dynamic phenomenon which has 
received a great deal of attention since the turn of the century. The concep- 
tual model of this problem is still relatively in its infancy. However, this 
much is known: that transition to turbulence can be caused by instability 

waves in the boundary layer. Modern airfoils are faced with a number of 
instabilities which can reduce the extent of laminar flow over the surface. 
Amplification of disturbances by Tol 1 mi en-Schl i chti ng , crossflow, or Gortler 
mechanisms over a certain threshold will cause premature transition. In order 
to maximize the extent of laminar flow, each, and subsequently any combination 
(the interaction problem) of these instabilities must be avoided or at the 
very least minimized. Of these, the TS and CF type disturbances are most 
commonly considered in modern airfoil design and analysis problems. 
Fortunately, computer codes designed to model these two disturbance types are 
available. Furthermore, correlation of these codes with experimental data is 
required to validate them as design tools. It is for these reasons that the 
TS and CF type disturbances are examined in the present analysis. 

Modern wings designed for transonic flight are generally swept. Figure 1 
shows a typical swept-back wing planform with its curved streamlines. Very 
near the surface, the boundary-layer streamline will have a curvature that is 
more pronounced than the potential flow streamline, which sets up a three- 
dimensional velocity field in the boundary layer. This three-dimensional 
field can be broken up into a streamwise component, parallel to the local 
potential streamline, and a crossflow component, perpendicular to the local 
potential streamline. Figure 2 shows the streamwise and crossflow components 
of a typical swept wing boundary-layer profile. Note that there is an 
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inflection point in the crossflow profile which is due to the imposed boundary 
conditions of no slip at the wall and the asymptotic approach to zero velocity 
at the boundary-layer edge. Rayliegh, as noted earlier, has shown that 
velocity profiles with inflection points are dynamically highly unstable even 
in the absence of viscosity. Thus, associated with the crossflow velocity 
component is an instability which is inviscid in nature and is due primarily 
to the existence of an inflection point in the velocity profile. This 
instability is referred to as crossflow (CF) instability. Velocity profiles 
without inflection points, for example the Blasius boundary layer, were 
initially thought to be inherently stable. Prandtl first noticed the mildly 
destabilizing effects of viscosity. Toll mien and Schlichting examined this 
type of profile in detail with the OS equation (which includes viscous 
effects) and also found viscosity to be destabilizing. Thus, the instability 
which arises from viscous effects was originally referred to as Toll mi en- 
Schlichting (TS) instability. However, for flow with non-zero pressure 
gradients (for instance our swept-wing case), inflection points can develop in 
the streamwise profile which is further destabilizing. An adverse or positive 
pressure gradient promotes highly destabilizing inflection points in the 
streamwise profile. Therefore the instability associated with the streamwise 
component can be due to either viscous effects or inflection points in the 
profile. 

Dagenhart, reference 13, discussed the fact that there are two wave angle 
orientations that maximize the disturbance amplification rates. These 
orientations correspond roughly with the streamwise and crossflow velocity 
directions. Squire has shown that for incompressible flow, the TS 
disturbances which amplify the most are those moving in the direction of the 
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free-stream flow. For this reason, the constraint of a wave angle of zero is 
applied to the TS stability calculations which follow. Gregory, Stuart and 
Walker (ref. 15) and Pfe.nn1n.ger (ref. 16) have shown by flow visualization 
techniques that CF vortices align approximately in the potential flow 
direction and are stationary. Furthermore, calculations by Srokowski and 
Orzag (ref. 17) indicate that these stationary disturbances are usually the 
most highly amplified. Therefore, a frequency of zero is assumed for the CF 
stability calculations which follow. 

Suction thins the boundary layer and modifies the profile shape. 17 
Figure 3 qualitatively illustrates the effect of suction on both the 
streamwise and crossflow profiles. In the case of the streamwise profile, 
suction draws the velocity profile down and tends to remove any inflection 
points. This greatly improves the stability of the streamwise flow. For the 
crossflow profile, however, suction is not as effective. As noted by Brown 
(ref. 18) and Mack (ref. 19), suction is more efficient at controlling TS 
instabilities than CF instabilities. The crossflow profile must have an 
inflection point because of the boundary conditions. Therefore, suction can 
only draw the inflection point down where viscous effects can influence and 
slow the disturbance growth rates. 

Computer Codes 

The incompressible stability characteristics of laminar compressible 
boundary layers on swept, tapered wings with suction are calculated by an Orr- 
Sommerfeld equation solver. This computational method can model both the TS 
and CF type disturbances within the boundary layer. In the analysis mode, two 
separate codes are required: one to calculate the boundary -layer velocity 

profiles, and one to calculate the stability characteristics. The boundary- 
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layer velocity profiles were obtained from a program written by Kaups and 
Cebeci (ref. 20) (referred to as the CEBECI code) while the stability 
characteristics were found using either the SALLY code of Srokowski and Orzag 
(ref. 17) or the MARIA code written by Dagenhart (ref. 13). 

The CEBECI code calculates the compressible three-dimensional velocity 
profiles of a laminar boundary layer on a swept and tapered wing with 
suction. An efficient two-point finite-difference method is used to solve the 
boundary-layer equations. Input to the program includes the streamwise 
pressure and suction distributions, free- stream flow conditions and wing 
geometry information. Upper and lower surface boundary layers are calculated 
separately. Program output which includes both the spanwise and streamwise 
components of the velocity distribution is saved in magnetic disk files in 
formats appropriate for both SALLY and MARIA codes. 

The SALLY code models the amplification of a travelling wave type 
disturbance from the initial point of amplification in a laminar boundary 
layer. The code utilizes Chebyshev polynomials to approximate solutions to 
the Orr-Sommerfeld equation. Input includes the boundary-layer velocity 
profiles from the CEBECI analysis along with the particular wavelength, 
frequency and/or wave angle of interest. Two of the three variables 
(wavelength, frequency and wave angle) may be specified while the third is 
part of the solution along with the amplification rate. Separate computer 
runs are required to model either the TS or CF type disturbances. For 
analysis of TS disturbances, the wave angle is fixed at zero (for 2-D waves) 
and the frequency range is specified. For analysis of CF disturbances, the 
frequency is fixed at zero (for standing waves) and the range of wavelengths 
is specified. Output is the natural logarithm of the growth of the 
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disturbance from the initial point of instability along the airfoil. It 
should be noted, however, that due to an inadequacy of the solution technique 
utilized within SALLY, the code is unable to calculate disturbance damping. 
This is due to the technique's ability to trace only eigenvalue solutions 
which predict amplification. 

The MARIA code approximates the amplification of the CF type disturbance 
only. A table look-up method based on stability charts generated by SALLY is 
used. This code has been shown (ref. 13) to adequately model the CF 
instability for preliminary and trade-off analyses and, if needed, shows the 
critical wavelengths for a more accurate SALLY run. The use of this code 
initially should reduce the overall computer time and costs. 
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LAMINAR FLOW CONTROL EXPERIMENT 


Design and Set-Up 

An experiment designed to evaluate the compatibility of laminar flow 
control and advanced supercritical airfoil concepts has been conducted 
recently by NASA in the Langley 8-Ft. Transonic Pressure Tunnel. This 
experiment examined the possibility of reducing friction drag through laminar- 
ization with local suction at simulated flight conditions which are 
representative of future transport designs. A large chord, swept 
supercritical airfoil with suction was tested over a wide range of flight 
conditions up to and within the transonic regime, which includes the design 
point case of M oo = 0.82 and R c = 20 million. A detailed discussion of the test 
design and setup was presented by Harvey and Pride (ref. 5), This LFC 
experimental concept was designed using the latest analytical methods. One of 
the major objectives of the experiment therefore was to validate these design 
techniques. The present results will not Include the experimental design 
point data due to a classification restriction on this information. However, 
the off-design test points will allow adequate examination of the disturbance 
growth problem. 

The experimental setup, as shown in figure 4, consists of the swept-wing 
model and the contoured liner, both of which were designed for tunnel 
conditions of = 0.82 and R c - 20 million and sweep angle A = 23°. The 
aerodynamic design considerations of the swept LFC wing were discussed in 
papers by Pfenninger, Reed and Dagenhart (ref. 21) and Allison and Dagenhart 
(ref. 22). The wing was designed to examine suction laminarization at 
conditions which simulate full-scale flight conditions. The design 
considerations of the contoured tunnel liner were discussed in a paper by 
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Newman, Anderson and Peterson (ref. 23). The contoured liner installed on all 
four walls of the test section is intended to simulate interference-free flow 
about an infinite-span, yawed wing at the design test condition. This is a 
point design for a highly sensitive transonic condition, therefore, contoured 
liner performance is expected to be degraded at off-design conditions. The 
effect of off-design testing will be discussed later. In order to approach 
free-flight disturbance levels within the test section, a combination of 
honeycomb, screens and choke-plates were installed in the tunnel. 

Measured Data 

Instrumentation for the LFC experiment includes pressure orifices, 
suction metering system, surface thin-film sensors for locating transition and 
acoustic transducers for measuring dynamic pressure data. These instruments 
provide the necessary information to correlate the n- factor method. The 
primary test surface of the LFC experiment is the upper surface, therefore, 
the present results are based on the upper surface data only. A sketch of the 
upper surface instrumentation is shown in figure 5. The pressure distribution 
along the densely-packed center-line row of orifices, plus the suction 
distribution are the basic input into the stability analysis in order to 
calculate the growth of disturbances within the boundary layer. This growth 
is then compared with the transition data in order to arrive at a transition 
n-factor. This is the basic method of correlation to be used in this paper. 

Figure 5 is a sketch illustrating both chordwise and spanwise rows of 
pressure orifices. This array of orifices provides the means for measuring 
the three-dimensional flow field about the upper surface of this airfoil. The 
flow field is expected to be non-uniform since the cases to be 
analyzed (M^ = 0.4, 0.6, 0.7) are at off-design Mach numbers. Shown in figure 
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6 for the center span are experimental chordwise pressure distributions 
(symbols) for all three Mach numbers against the theoretical design point 
pressure distribution (dashed line). This figure reveals that the 
experimental pressure distribution differs considerably from design for all 
three Mach numbers in the following regions: (1) 0.2 < x/c < 0.20 (a steep 

adverse pressure gradient), and (2) 0.20 < x/c < 0.80 (a slightly favorable 
pressure gradient). In the region where there is an adverse gradient, 
disturbance amplification is expected, while damping is likely for a region 
with a favorable gradient. Since the pressure distribution is scaled by the 
free-stream dynamic pressure, the three Cp plots appear similar. Close 
inspection, however, shows that the pressure peak rises slightly with 
increasing Mach number. Also, as can be seen in figure 6, these three Mach 
numbers are all subcritical. Inspection of the spanwise pressure 
distributions shown in figure 7 shows that a gradient exists along lines of 
designed constant pressure. These design isobars coincide with the suction 
slots running along lines of constant percentage of chord. The spanwise 
gradient effect is more pronounced at the leading and trailing edges. The 
fact that the pressure distribution changes with spanwise location suggests 
that as the suction levels are changed the transition pattern will move non- 
uni formly. This phenomenon will be seen in the transition data to be 
presented later. This spanwise gradient is an indication that the liner 
performance is degraded at these off-design Mach numbers. 

Suction air is removed through closely spaced fine slots extending span- 
wise along constant percentage of chord. Details of the suction system design 
are given in reference 5. Evaluation of the suction rates required to avoid 
premature transition on the swept LFC airfoil was based on incompressible 
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stability analysis at the design condition. Since non-ideal pressure distri- 
butions are being analyzed, the optimum suction distribution required for 
laminarization differs from the design suction. This is shown in figure 8, 
where two theoretical suction distributions, one for the design case and one 
for an off-design case, are plotted with two experimental suction distri- 
butions. It can be seen that the experimental suction levels are quite a bit 
higher than the theoretical levels. The experimental suction distributions 
required to maintain fun-chord laminar flow are higher because the Mach 
numbers are off-design. 

Suction levels required to maintain laminar flow are greatly influenced 
by Reynolds number. A parameter, which is referred to as Reference Suction 
Level (RSL), was defined in order to remove this Reynolds number dependence. 
The RSL parameter is defined as 


where 


RSL = C p /C p 


REF 


C Q = /J C p (x/c)d(x/c) 


and 





'DES 


C f . (x/c)d(x/c). 
W DES 


A detailed discussion of the development of this parameter is in the 
Appendix. The RSL coefficient is the actual integrated suction level scaled 
by a reference integrated suction level which is based on Reynolds number. 
This essentially allows suction parameter comparisons based on removal of 
equal fractions of the boundary layer. 

Throughout the experiment, the suction distribution maintained the same 
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basic shape as the design distribution-only the overall levels were changed. 
This allows the RSL parameter to be an indication of the effect of overall 
suction on boundary-layer instabilities. Experimental data were taken by 
setting the desired Mach and Reynolds number combination in the test section 
with overall suction at its maximum level and then examining the transition 
pattern as the overall level was reduced. This provided several test cases 
for a range of suction levels at each Mach and Reynolds number. This is 
referred to as a suction excursion. As will be discussed later, the suction 
excursion is essential in order to arrive at transition n-f actors. Thus, a 
RSL parameter variation, which represents a suction excursion, will be the 
primary method used to examine the dependence of boundary-layer stability on 
suction levels. 

Hot-film sensors deposited on one end of a tiny quartz rod were mounted 
flush with the surface at locations shown in figure 5. The sensors were 
controlled by a constant temperature anemometer at 80°C above local recovery 
temperature. The installation technique for the hot-film sensors on this 
model is given in reference 5. The hot-film sensors were used to analyze the 
local boundary-layer state. There were a total of 25 sensors distributed on 
the upper surface: 9 along the mid-chord line, with the remainder positioned 

in a dense array in the aft pressure recovery region. This arrangement was 
chosen to provide detailed coverage of the changing transition pattern at the 
design condition, which was expected to move in a two-dimensional fashion. 

For off-design conditions, however, the transition pattern moved non- 
uni formly. Therefore the distribution of sensors shown was not as thorough as 
was desirable. The effect of this will be discussed later. Examples of some 
voltage traces, their corresponding boundary- layer states and their locations 
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on the upper surface are shown in figure 9. The transition pattern shown is 
typical of the data analyzed. The amplitude of the fluctuating value of the 
sensor voltages indicates the relative condition of the boundary layer. The 
low amplitude dynamic signal indicates the lower shear stress that is 
representative of a laminar boundary layer and is typically found as shown in 
the figure. The higher amplitude dynamic signal indicates the higher shear 
stress that represents a turbulent boundary layer. The transitional type 
boundary layer has both the laminar and turbulent type signals and it is the 
relative amount of each, or intermittency factor (IF), which indicates the 
fraction of time that the flow is turbulent. For example, a signal that shows 
turbulent flow of 80 percent of the time has an IF of 8. 

A major component of a complete stability analysis, as noted in recent 
studies by Spangler and Wells (ref. 24) and Mack (ref. 25), is an under- 
standing of the disturbance environment influencing transition. Mack notes in 
reference 25 that any empirical transition criteria can only be valid for a 
specific disturbance environment. Therefore, measurements of the dynamic 
pressure data were taken in order to assess the specific disturbance spectral 
levels that were present in the test section for different test conditions. 
Acoustic transducers, which are Hemholtz-cavity mounted within cylindrical 
elliptic-nose probes, were used to measure the fluctuating pressures. The 
probes were mounted on a strut located upstream of the wing leading edge in 
order to measure the free-stream disturbance levels impinging on the model. 
Dynamic signals were recorded on a tape recorder, and were later analyzed 
using a high resolution digital spectrum analyzer. The transducer signals 
were corrected for electronic noise using a technique outlined in reference 
26. A sample energy spectrum is shown in figure 10 with electronic noise 
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subtracted out of the raw data and then converted to obtain the corrected 
disturbance spectrum. 
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STABILITY ANALYSIS 


The upper surface LFC experimental boundary layer was analyzed for TS and 
CF disturbances. Only the off-design results are presented. The TS analysis 
was completed using the SALLY code, while the CF analysis was completed using 
the MARIA code. A total of 90 separate data cases were analyzed in which the 
3 main parameters (RSL, M^ and R c ) were varied. To show representative n- 
f actor calculations, one set of data cases was chosen in which the suction 
parameter was varied while the Mach and Reynolds numbers were held constant. 
This set was chosen to illustrate how the transitional n-factors were inferred 
from the analyzed data. Summary plots will follow in which the results of all 
the cases are presented. Attached laminar flow up to the experimentally 
determined transition location was computed for all cases considered. 

Suction Excursion 

As stated before, the transition pattern moves non-uni formly spanwise as 
the overall suction level is changed for the off-design Mach number cases 
being analyzed. The reason for this is that the tunnel liner is mis-matched 
for its intended purpose at these Mach numbers. The liner was intended to 
provide infinite swept flow (the chordwise pressure distribution should remain 
constant along the span) at the design condition. Since the off-design 
conditions are being analyzed, the liner does not correctly divert the flow in 
an infinite swept fashion. Therefore, the transition pattern moves three- 
dimensionally. Figure 11 shows that as the suction levels are reduced for 
constant M and R , as indicated by a reduction in RSL, the transition pattern 
sweeps down towards the center-span row of orifices. At some point, when the 
suction level is reduced far enough, the line corresponding to the beginning 
of transition jumps forward at mid-span. Along the center row of orifices. 
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the transition pattern almost instantaneously changes from 100 percent chord 
laminar to nearly 100 percent chord turbulent. Because of this phenomenon, a 
different method of correlating n- factors with transition was used for this 
analysis. In previous methods, the transition n-f actor was found by 
calculating an n-f actor corresponding to the experimentally determined 
transition point. For the present method, transition n-factors were found by 
locating the maximum calculated n- factor in the nose region at the moment 
transition moves forward of the initial thin-film gage. Thus, it takes a 
suction excursion of points to estimate the correct n- factor corresponding to 
transition at each Mach and Reynolds number combination. 

Sample Cases 

Inferred transitional n-factors were arrived at by examining suction 
excursions at each Mach and Reynolds number combination. To illustrate the 
basic correlation method, a representative excursion will be discussed in some 
detail. The cases that follow represent a suction excursion at M^ = 0.6 and 
R c = 10 x 10 6 . Figures 12 thru 15 present the results of the SALLY and MARIA 
calculations for each of the 4 cases. Shown on the plots are the calculated 
logarithmic amplification ratios for a range of frequencies or wavelengths and 
the pressure and suction distributions along the chord. For the curves 
showing the T$ instability, the dotted lines represent the area of neutral 
stability or damping and the symbols correspond to the frequencies listed in 
the table. For the CF curves, only the most critical wavelength and 
corresponding maximum n-factors are listed. 

Case 1 

The TS and CF stability calculations for a suction parameter RSL = 1.8 
are shown in figure 12. The corresponding transition location along the 
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center station, as shown in figure 11, is (x/c) tr > 95 percent. Figure 12a 
shows that for TS growth, the most amplified frequency is f = 14 kHz which 
reaches a logarithmic amplification ratio of n = 7.413 in the nose region, 
followed by a mid-chord stable or damped region and a subsequent growth near 
the trailing edge. The growth that is calculated near the tail occurs past 
the last suction slot and is beyond all thin-film instrumentation. 

Furthermore, since the mid-chord region, in which no amplification is 
predicted, is large and since SALLY is unable to predict damping, the actual 
disturbance amplitudes near the trailing edge should be much lower. For these 
reasons the growth near the tail is ignored. As stated in figure 12b, CF 
disturbances only experience an integrated logarithmic amplification of 
n = 0.72 in the nose region for the critical wavelength x/c = 0.0006, followed 
by a stable region to the tail. Obviously, TS amplification is the dominant 
instability. Furthermore, a TS growth of n = 7.4 at the nose did not force 
transition to occur anywhere along the chord. 

Case 2 

The stability calculations for a suction parameter RSL = 1.37, with a 
corresponding full -chord laminar flow, are shown in figure 13. For TS growth, 
figure 13a shows that for the most amplified frequency f = 14 kHz, a 
logarithmic amplification ratio of n - 9.0 is reached in the nose region, 
followed by a mid-chord damped zone and then a subsequent growth near the 
tail. Again, the growth in the tail region is overpredicted and is ignored. 

As depicted in figure 13b, the maximum amplification ratio for CF instability 
is n = 0.78 for a critical wavelength X/c = 0.0006 and this occurs in the nose 
region. TS amplification is still the dominant instability and a growth of 
n = 9.0 in the nose region still did not produce early transition. 
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Case 3 

Figure 14 shows the stability calculations for a suction parameter RSL = 
1.25. As can be seen in figure 11, the transition location along the center 
span has moved forward slightly to (X/C) tr « 85 percent. Figure 14a shows 
that for TS growth, the most amplified frequency is f = 11 kHz which reaches a 
logarithmic amplification ratio of n = 9.84 in the nose region, followed by a 
mid-chord stable region all the way back to transition. For CF instability, a 
maximum amplification ratio of n - 0.76 is reached at the nose for a critical 
wavelength of x/c = 0.0006 as stated in figure 14b. The CF growth is still 
relatively insignificant. However, a TS growth of n = 9.84 in the nose region 
corresponds to a 15 percent forward shift in the transition location. 

Case 4 

The TS and CF stability calculations for a suction parameter RSL = 1.15 
are shown in figure 15. The transition location corresponding to this case, 
as shown in figure 11, has moved forward of the initial thin-film gage. 
Therefore transition has occurred fairly close to the leading edge. Figure 
15a shows that for TS growth near the leading edge, the most amplified 
frequency is f = 12 kHz which reaches a logarithmic amplification of n = 10.52 
before becoming stable. Figure 15b depicts that for CF growth near the 
leading edge, the most amplified wavelength is x/c = 0.0006 which reaches a 
logarithmic amplification ratio of n = 0.77. TS amplification remains the 
dominant instability. More importantly, a TS growth near the nose of 
n = 10.52 corresponded with transition occurring near the leading edge. 

Data Summary 

Stability analyses of off-design LFC experimental data, over a range of 
Mach numbers, Reynolds numbers, and suction levels, was completed and the 
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results are tabulated in table 1. TS instability was found to dominate the 
process leading to transition while CF growth was small enough such that 
interaction of these two disturbances seems unlikely. A summary of the TS 
amplification for all cases considered is presented in figures 16, 17 and 
18. Each figure corresponds to one Mach number in which the calculated TS 
logarithmic amplification ratios (ny$) are plotted against Reynolds number 
(R c ) for a number of suction levels (RSL) . The symbols represent the 
condition of the flow - open symbols for full -chord laminar flow, closed 
symbols for full -chord turbulent and half-open symbols for transition occuring 
along the chord between the most-forward and most-rearward thin-film gages. 
These plots are used to indicate at what level of TS amplification transition 
is likely to occur. The dashed lines that are sketched in indicate the n- 
f actor limits corresponding to transition. The symbols with flags signify a 
phenomenon, referred to as oversuction, that was observed in the transition 
data. This phenomenon shall be discussed in more detail later. The three 
data sets corresponding to each Mach number were not as complete as was 
desirable. At the lowest Mach number, the range of suction test cases 
available for analysis did not provide a sufficient band of n-factors at 
transition. Therefore, the figures will be discussed in descending order of 
completeness. 

A phenomenon was observed in the transition data for the cases with the 
combination of high suction levels, and high Mach and Reynolds numbers. This 
phenomenon amounted to a reduction in the laminar zone at the highest suction 
levels which disappeared as the suction levels were reduced. Shown in figure 
19 is an example of the transition data which demonstrates this phenomenon. 

The figure depicts transition patterns for cases with a range in suction level 
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(RSL) at = 0.6 and R c = 18 million. The lines shown on the upper surface 
plot represent the beginning of transition for each case as revealed by the 
hot-film sensors. At the highest suction levels (RSL = 2.56) transition 

occurs as far back as 100 percent on either side of mid-chord, but along the 

mid-chord line transition is at about 85 percent chord. As RSL is reduced to 

2.16, the transition line moves towards the tail and is at ( x/c ) « 90 

percent along the mid-chord. For RSL = 2.01 transition at the mid-chord line 
has moved back beyond the last hot- film sensor, but the transition line has 
moved forward for the outboard region at the top of the figure. As the 
suction is reduced even further, transition again jumps forward as discussed 
earlier and as shown in figure 19, It seems that this loss of laminar flow at 
high suction levels may be due to sensitivity to surface roughness. 

Inspection of figure 5 which shows the instrumentation (pressure orifices, 
hot-film and acoustic sensors) suggests that the maximum roughness due to 
instrumentation occurs along the centerline of the model. Too much suction 
will excessively thin the boundary layer, such that transition is governed by 
surface roughness and not by disturbance amplification. These cases that 
demonstrate such behavior are referred to as oversuction cases and are shown 
in figures 16, 17, and 18 as half-open symbols with a flag. 

The most complete data set is at = 0.6. The data shown in figure 17 
clearly define the transitional limits. Suction levels which allow TS 
amplification factors greater than n » 10 at the lower Reynolds numbers 
to n * 11 at the higher Reynolds numbers are insufficient to maintain laminar 
flow. On the other hand, suction levels which allow up to n « 9 at the lower 
Reynolds to n * 10 at the higher Reynolds numbers should maintain full -chord 
laminar flow. The effect of surface roughness can be seen in the figure. 
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where at the higher Mach and Reynolds numbers, suction must be reduced to 
attain full -chord laminar flow. To summarize at = 0.6, logarithmic 
amplification of TS type disturbances with n > 10 leads to transition. 

The next best data set is at M^ = 0.7. Figure 18 shows that the data 
available indicate the allowable n-factor range at transition. A TS 
logarithmic amplification of n « 10 is allowable without reduction in full- 
chord laminar flow. However, TS amplification ratios of n > 10.5 for the low 
Reynolds numbers to n > 11 for the high Reynolds numbers indicates that 
transition will occur far forward on the model (essentially full-chord 
turbulent flow). The correlation trends for this set of data is in close 
agreement with the data at = 0.6. However, the effect of surface 
roughness appears to be more critical at M^ = 0.7. This can be seen in 
figure 18 where the flagged data cases are encroaching on the transition 
limit, to the point where full-chord laminar flow may be difficult to attain 
at any suction level for the higher Reynolds numbers. Excessive boundary- 
layer thinning may be a problem which warrants further study at the higher 
Mach and Reynolds numbers. 

The last data set is for = 0.4 and is shown in figure 16. Because, at 
this Mach number, it was difficult to get transition along the center span 
even at the minimum attainable suction levels, the data seem inconclusive. 
However, based on the results at the other two Mach numbers and some 
engineering judgment, transition limit lines can be sketched in as shown. If 
this speculation is correct, then a TS logarithmic amplification of n « 8.5 at 
the low R c to n « 9 at the high R c would correlate with transition. This 
trend is in agreement with the higher Mach numbers, even though the levels are 
a bit lower. For this Mach number, surface roughness effects appear in the 
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transition data only at the highest Reynolds number, R c = 20 million. 

Discussion of Results 

The off-design conditions induced .spanwi.se pressure gradients and large 
minimum pressure peaks near the nose of the airfoil. This produced spanwise 
variations in boundary-! ay er development, disturbance amplification near the 
leading edge and a large stable TS zone. Three-dimensional pressure and 
transition patterns developed which caused the abrupt jump of transition to 
the leading edge as the suction level was reduced. All these factors 
contributed in increasing the difficulty of correlating n-factors at 
transition. However, despite the non-ideal conditions, the calculated n- 
f actors at transition are in reasonable agreement with previous studies. 

The TS disturbances were amplified only over a small region near the 
leading edge due to a combination of factors. First, sharp minimum pressure 
peaks were generated near the leading edge because the Mach numbers considered 
were less than the design Mach number. This creates a negative pressure 
gradient, which is greatly destabilizing for TS disturbances, and therefore, a 
region where, in the absence of suction, strong disturbance growths could be 
expected. As stated before, TS disturbances are greatly stabilized by 
suction. Secondly, the suction distribution and, more importantly, the 
suction system (slot geometry, spacing and locations) were designed for a 
free-stream Mach number of 0.82 and Reynolds number of 20 million. The 
initial suction slot is located at approximately 2.5 percent chord which is 
slightly aft of the pressure peak. In general, TS disturbances have begun to 
amplify before the initial suction slot. Suction should, ideally, be applied 
just prior to and over a region of disturbance growth. Thus, in the sense 
that amplification is allowed to take place without suction being applied, 
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this present suction system is not ideal for handling the off-design Mach 
numbers. The analysis also shows that, once suction is applied, TS 
disturbances are damped very quickly and remain stable past the last suction 
slot. The entire amplified band of frequencies is stabilized, for all suction 
levels, roughly 2.5 percent chord after the initial suction slot. This 
implies that if suction were available closer to the leading edge and just 
prior to the pressure peak, then the calculated growth rates for these off- 
design cases would be less severe. Thus, full -chord laminar flow could be 
maintained for both lower values of suction and higher chord Reynolds numbers 
at each Mach number. 

The calculated logarithmic amplification ratios at transition change with 
Mach number, as depicted in figure 20. Both TS and CF n-factors at transition 
are shown for the range of Reynolds numbers analyzed. Based on the figure, 
one can see that TS amplification is the dominant instability. CF 
amplification is relatively small and does not appear to be a problem until 
the Mach number approaches the design value. Curiously, the TS transitional 
n-factors are increasing slightly with Mach number. Examination of the noise 
levels, in terms of measured pressure fluctuations (in percent of local free- 
stream conditions), for the 8-Foot Transonic Pressure Tunnel shown in figure 
21 (taken from ref. 27), shows that acoustic noise levels increase with Mach 
number. One might expect that based on this, the transition n-factors would 
decrease with Mach number. Since the opposite happens, perhaps the overall 
noise levels are relatively unimportant. Figure 22 shows the disturbance 
frequency spectra for the Mach numbers in question. These sample curves, 
taken from the acoustic data measured upstream, show that the overall noise 
level, represented by the area under each curve, increases with Mach number. 
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However, in the range of frequencies critical for TS amplification (10 < f < 

20 kHz), the curve for = 0.4 has peak fluctuations higher than that at the 
other two Mach numbers. Therefore, it appears that the disturbance amplitudes 
within the critical band of amplified frequencies is more important than the 
overall sound pressure level at influencing transition. 
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CONCLUSIONS AND RECOMMENDATIONS 

Stability analysis of existing experimental data available from the NASA 
Langley 8-Foot LFC experiment was completed. As a result of this analysis a 
number of conclusions and recommendations can be made: 

1. For the off-design cases analyzed, incompressible linear stability 
theory correlates with a n-f actor for TS waves of approximately 10 at 
transition. 

2. Tollmien-Schlichting amplification was the dominant instability 
influencing transition for all the cases analyzed while crossflow 
amplification was negligible. 

3. The off-design conditions provided a unique growth pattern in which 
disturbance amplification was limited to a narrow region near the nose 
and a stable region dominated the remaining 95 percent chord. 

4. The Tollmien-Schlichting n-factors at transition increased slightly as 
the Mach number increased, which may be due to the changing 
disturbance spectra within the narrow band of amplified frequencies. 

5. The good correlation between n-factors and transition over the broad 
range of non-ideal test conditions indicates that this method is a 
reliable design tool for advanced LFC airfoils in subsonic flight. 
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6. Future LFC wing designs should consider incorporating suction as far 
forward as possible with the capacity to significantly increase the 
suction for off-design conditions such as those studied here. 

7. Results of this off-design analysis indicate that transition due to 
surface roughness, and not disturbance amplification, is a problem 
that increases as the Mach and Reynolds number increase. This 
oversuction effect could be significant at the design case. Further 
studies of this phenomenon at the design case is recommended. 

8. Future consideration is recommended for examination of compressibility 
effects on disturbance amplification. 
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APPENDIX 


Development of RSL Parameter 
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The laminar boundary-layer thickness is known to vary with Reynolds number 
as (ref. 14) 

6/c~l /nC Q . (1) 

In order to remove the same fraction of the boundary layer for various test 
Reynolds numbers, the suction coefficient should also vary as the boundary- 
layer thickness varies 

Cq ~ 5/c ~ 1 / /" R” . (2) 

Therefore, for the same fraction of the boundary layer removed for a test case 
compared to a design case we have 



The design integrated suction coefficient is defined as 

C n = ji C r (x/c)d(x/c) . (4) 

W DES u %S 

The object is to remove the Reynolds number dependence from the suction 
parameter. Therefore the integrated suction coefficient for the case of 
interest shall be compared to a reference integrated suction value that is 
scaled by the Reynolds number. Thus the Reference Suction Level (RSL) 
parameter is defined as 
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Rsl = e Q / t Q 


(5) 


REF 

where the integrated suction coefficient of interest is 

Cq = jJ Cq(x/c)d(x/c) (6) 

and the reference integrated suction level is the design suction level scaled 
by the Reynolds number difference 




/ jr —rw c 

C DES c 


(7) 


The design suction distribution and overall level is shown in figure 8. Based 
on the integrated suction value for the design case, equation (7) was used to 
generate the reference levels as shown in the table 


R 


c 



10 million 

-0.00049 

12 

-0,00045 

14 

-0.00042 

16 

-0.00039 

18 

-0.00037 

20 

-0.00035 


These reference values were used in conjunction with equation (5) to find the 
RSL for each case. 
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Table 1.- Tabulated results of incompressible stability analysis. 
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Table 1.- Continued 
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Table 1.- Concluded. 
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Figure 1.- Flow conditions on typical sweptback wing* 
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Figure 2.- Swept-wing boundary-layer profile. 
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SUCTION FLOW 
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Figure 4. -Schematic of the tunnel and model for LFC experiment in the 
NASA Langley 8-ft Transonic Pressure Tunnel. 
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Figure 5.- Schematic of instrumentation locations on upper surface of 
LFC model . 
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x/c 

Figure 8.- Comparisons of theoretical and experimental suction distributions. 



Turbulent: 
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Figure 9.- Example hot-film signals and corresponding locations on upper 
surface of LFC model. 




Peak disturbances 



Figure 10.- Sample free-stream acoustic data for LFC test 
Moo =.4 and R =10 million. 
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Figure 11.- Effect of suction excursion on transition pattern 
Mo,, =. 6 and R-=10 million. 
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Figure 12(b).- Incompressible CF calculations for case 
R_ = 10 million and RSL = 1.84. 
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Figure 13(a).- Incompressible TS calc 
R r = 10 million and RS 
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Figure 13(b).- Incompressible CF calculations for case 
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iressible TS calculations for case 
0 million and RSL = 1.15. 
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Figure 15(b).- Incompressible CF calculations for case 
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Chord Reynolds Number, 
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Chord Reynolds Numbers, 
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Chord Reynolds Number, 
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Figure 19.- Effect of oversuction on transition. 
Mo =.6 R c =18 million 


68 



A = 23? upper surface 
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Figure 20.- Calculated incompressible n-factors at transition 
for swept LFC airfoil. 
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Figure 21.- Measured pressure fluctuations in the 8-Ft. TPT (ref. 27). 
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14x10 




Figure 22.- Sample noise spectra from upstream acoustic transducers. 
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